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We complete our analysis of Nf=2 QCD based on the lattice staggered fermion formulation. Using 
a series of Monte Carlo simulations at fixed miL\ h one is able to test the universality class with 
given critical exponent y/,. This strategy has been used to test the 0(4) universality class and it 
has been presented at the previous Lattice conferences. No agreement was found with simulations 
in the mass range am q = [0.01335,0.15] using lattices with L s = 16 up to 32 and L t = 4. With 
the same strategy, we now investigate the possibility of a first order transition using a new set of 
Monte Carlo data corresponding to yt, = 3 in the same mass and volume range as the one used for 
0(4). A substantial agreement is observed both in the specific heat scaling and in the scaling of 
the chiral condensate, while the chiral susceptibilities still presents visible deviation from scaling 
in the mass range explored. 
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1. Motivation 

The order of the chiral transition of two-flavor QCD is still debated. Its relevance is both 
phenomenological (heavy ion experiments) and theoretical: is deconfmement an order-disorder 
transition associated with a change of symmetry? 

With massless quarks the chiral/deconfinement transition for N/=2 is expected to be second 
order in the 0(4) universality class or first order. Two completely different scenarios correspond 
to those two possibilities. Since, contrary to the first order case, second order phase transitions are 
unstable against the explicit breaking of the underlying symmetry, in the second case one would 
have a crossover instead of a real phase transition for small but non-zero quark masses. That would 
imply the possibility of going continuously from a confined to a deconfined state of matter, in 
contrast with the idea of confinement being an absolute property of strongly interacting matter at 
zero temperature and of deconfmement being an order-disorder transition associated to a change of 
symmetry [Jp. A second consequence would be the presence of a crossover line, at finite mass, in 
the region of high temperature and small baryon chemical potential }Xb of the QCD phase diagram, 
thus implying a critical point [§] to connect with the first order line which is supposed to exist at 
low temperatures and large chemical potentials. No such point is expected to exist if the transition 
at = is first order. No critical point has been found up to now in experiments with heavy ions, 
but the question is still open I@» □» §1- 



2. The method 

To establish the order of a phase transition using Monte Carlo simulations is not an easy task. 
Since all quantities are analytical in a finite volume, one must study the finite size scaling (FSS) 
to look for developing singularities. Here the FSS analysis is made intricate by the fact that, since 
simulations at zero quark mass are not feasible, one has to consider both the finite spatial size L s 
and the finite quark mass mi = am q . The FSS ansatz for the singular part of the free energy density 
F s around the chiral critical point is: 

■l/v 



F s (z,m L ,L s ) = L~ F s [%L l s '\m L Ll^ , (2.1) 

where % is the reduced temperature T = (1 — T/T c ), V and yu are critical indexes. From the above 
equation we can derive the FSS of the specific heat: 

C v - Co ~ L? /v > c (tZ,! /v , mJJ^ , (2.2) 

and of the susceptibility % m of the chiral condensate: 

Xm-Xo^ Lj /V <l> x (xLl /v ,m L LyA . (2.3) 



In Ref. ^ we have approached the problem varying and L s in such a way that the scaling 
variable rtiiLl 1 ' was fixed, thus reducing the FSS Eqs. (2.2) and (2.3), to one variable. Since the 



critical exponent yj, must be fixed, one can only test a given universality class, which in Ref. [g] 
was chosen to be 0(4). In that case no sign of scaling with the 0(4) critical exponents was found. 
In the present work we repeat the analysis fixing v/, = 3 as expected for a first order transition. 
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3. Simulation details 

As in Ref. [^], the lattice discretization used is the standard staggered action. The algorithm 
used in the new simulations for the test of first order scaling, is the RHMC. All lattices have 
L t = 4. New Monte Carlo data have been generated to match the L s = 32, mi = 0.01335 lattice 
used in Ref. [^], thus fixing the parameter miJ? s . New datasets for lattices with spatial extension 
L s = 16,20,24 have been produced. For each of them ten different values of j8 spanning the entire 
critical region were simulated with a total statistics of about 90k thermalized trajectories collected 
for each of the lattices. The data collected is analized using a multi-histogram reweighting to 
extrapolate the observables to intermediate couplings. The backgrounds Co and % appearing in 



Eqs. ( |2.2| ) and ( [2.3D are determined as in Ref. [|6p, i.e. Co from a linear fit of the tails of the measured 
plaquette susceptibilities, and Xo is fixed by a fit of the maxima of the chiral susceptibility in the 
critical region. 



4. Test of scaling 

In Fig. [T] the scaling expected for a first order transition for the specific heat Cy and of the 
chiral susceptibility Xm i s shown. The picture also includes the corrisponding figures for 0(4) 
from Ref. [^]. For the first order case a reasonable scaling is observed for the specific heat Cy. 
As remarked in Ref. [Qj, Cy is independent of any prejudice on the symmetry and on the order 
parameter. Violations of the scaling Eq. ( [2J| ) are observed for Xm at larger values of the masses. In 
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Figure 2: Comparison of first order and (9(4) scaling of the subtracted chiral condensate. 



fact Eq. (23) is expected to be valid for the susceptibility of the order parameter. At large masses 
chiral symmetry is badly broken and is possibly not a good order parameter. In our data, 

both of the present paper and of Ref. [[J, Eq. ( |2.3[ ) seems to be violated for ml > 0.05. 

The scaling of the chiral condensate (magnetic equation of state) can also be checked. The 
comparison between first order and second order 0(4) is shown in Fig. |2[ The first order scaling 
is very good. As in Ref. [^], the scaling of the chiral condensate is obtained after a subtraction 
proportional to m L as m L — > 0. No sign of scaling is observed for the second order case. 



5. Conclusions 



The direct test of first order scaling shows a good scaling for the specific heat and for the chiral 
susceptibility at masses mi < 0.05. At larger masses scaling of the chiral susceptibility is broken, 
presumably because of strong breaking of the chiral symmetry. By comparison with the similar 
test for 0(4), the first order scaling is clearly preferred over the second order. We can say that 0(4) 
is inconsistent with the lattice data while a first order transition describes well the observations. 
Possible effects due to the discretization must be taken into account. We believe that ultraviolect 
effects should be irrelevant with respect to the large volume behavior. However the use of finer 
lattices and new simulation algorithms to approach the chiral limit will possibly clarify this issue. 
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